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Abstract 

Critical properties of hulls of directed spiral percolation (DSP) clusters are studied on the square 
and triangular lattices in two dimensions {2D). The hull fractal dimension (dn) and some of the 
critical exponents associated with different moments of the hull size distribution function of the 
anisotropic DSP clusters are reported here. The values of dn and other critical exponents are found 
the same within error bars on both the lattices. The universality of the hull's critical exponents then 
holds true between the square and triangular lattices in 2D unlike the cluster's critical exponents 
which exhibit a breakdown of universality. The hull fractal dimension {dn ~ 1.46) is also found 
close to 4/3 and away from 7/4, that of ordinary percolation cluster hull. A new conjecture is 
proposed for the hull fractal dimension (dn) in terms of two connectivity length exponents (i^y 
& z/_l) of the anisotropic clusters generated here. The values of dn and other critical exponents 
obtained here are very close to that of the spiral percolation cluster hull. The hull properties of 
the DSP clusters are then mostly determined by the rotational constraint and almost independent 
on the directional constraint present in the model. 



I. INTRODUCTION 



Directed spiral percolation (DSP), a new site percolation model, is recently introduced 
by Santra and Sinha . The DSP model is constructed imposing both directional and 
rotational constraints on the ordinary percolation model 0]. The directional constraint is in 
a fixed direction in space and the empty sites in that direction are accessible to occupation. 
Due to the rotational constraint the sites in the forward direction or in a rotational direction, 
say clockwise, are accessible to occupation. The direction of the rotational constraint is not 
fixed in space and it depends on the direction from which the present site is occupied. 
The cluster properties of the DSP model has already been studied on the square jl| and 
triangularis lattices in two dimensions {2D) at their respective percolation thresholds. The 
DSP clusters were found highly rarefied and anisotropic with chiral dangling ends. It was 
observed that the values of the critical exponents obtained on both the lattices are different 
from that of other percolation models like ordinary percolation (OP)^], directed percolation 
(DP)!^], and spiral percolation (SP)j^. The clusters on the triangular lattice were found 
more compact and less anisotropic than the clusters on the square lattice due to higher 
number of connectivity on the triangular lattice. It was found that the DSP model not 
only belongs to a new universality class than other percolation models but also exhibits a 
breakdown of universality in the cluster properties between the square and triangular lattices 
in 2L)0. 

The hull of a cluster is a continuous path of occupied sites at the external boundary 
of the cluster. The hull has important significance of its own apart from being a part of 
the percolation cluster. It has connection with the physical problems like nucleation and 
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growthp, surface reaction[7|, diffusion[8| etc. The hulls also exhibit scaling behaviour, dif- 
: erent from that of clusters, characterized by critical exponents at the percolation threshold 

in 

91. Ilfl|. Since the clusters generated in the DSP model are anisotropic and chiral in nature, 
it is interesting to study the critical behaviour of these anisotropic hulls with chiral dangling 
ends at p = Pc- 

In this paper, the critical properties of hulls of anisotropic DSP clusters are studied on 
both the square and triangular lattices. Results of anisotropic hull are not reported before 
and this is the first study of anisotropic-chiral hull properties. A comparison is made between 
the results obtained on the two lattices as well as in different models, like OP and SP. 



II. CLUSTER GENERATION AND HULL EXTRACTION 

A single cluster growth Monte Carlo (MC) algorithm has already been developed in Ref. P, 
0] to generate clusters under the presence of both directional and rotational constraints 
following the original algorithm of Leath[ll]. A brief description of the model is given 
here. In this algorithm, the central site of the lattice is occupied with unit probability. 
All the nearest neighbours of the central site can be occupied with equal probability p in 
the first time step. As soon as a site is occupied, the direction from which it is occupied 
is assigned to it. Selection of empty nearest neighbours in the next MC time steps is 
illustrated for the square lattice in Figure ^a) and triangular lattice in Figure Wib) . The 
descriptions given below are valid for both Figure Wict) and Wib) . Two long arrows from left 
to right represent the directional constraint. The presence of the rotational constraint is 
shown by the encircled dots. The black circles represent the occupied sites and the open 
circles represent the empty sites. The direction from which the central site is occupied 
is represented by a short thick arrow. The eligible empty site for occupation due to the 
directional constraint is indicated by the dotted arrow and the thin solid arrows indicate the 
eligible empty sites for occupation due to the rotational constraint. There are two empty 
sites due to the rotational constraint eligible for occupation at any MC step on the square 
lattice, as shown in Figure ^a) , whereas on the triangular lattice there are three such sites 
available, as shown in Figure HJ^fe) . After selecting the eligible sites for occupation, they are 
occupied with probability p. The coordinate of an occupied site in a cluster is denoted by 
{x,y). Periodic boundary conditions are applied in both directions and the coordinates of 
the occupied sites are adjusted accordingly whenever the boundary is crossed. At each time 
step the span of the cluster in the x and y directions = Xmax — Xmin and Ly = ymax —ymin 
are determined. If or Ly > L, the system size, then the cluster is considered to be a 
spanning cluster. The critical percolation probability pc is defined as below which there is 
no spanning cluster and at p = pc a spanning cluster appears for the first time in the system. 

As soon as a cluster is generated, the hull of the cluster is determined by the method 
of Ziff et al [12j. The hull of a cluster is defined to be the continuous path of occupied 
sites at the external boundary of the cluster. To determine the hull, a pair of occupied and 
empty nearest neighbour sites on the external boundary is chosen. An arrow is drawn from 
the empty site to the occupied site to define a direction and one moves to the occupied 



site. Facing to the direction of the arrow, a search is made starting from the left for an 
occupied nearest neighbour. As soon as an occupied site is encountered an arrow is drawn 
from the present site to the new occupied site and one moves to the new occupied site. 
The process is continued until the path passes the starting point in the same direction as 
it had first been passed. All the occupied sites encountered during this walk are listed in 
an array called as hull of the cluster. Typical hulls extracted from the clusters generated at 
the percolation thresholds on the square and triangular lattices are shown in Figure |2l and 
Figure El respectively. It could be seen from the hull itself that the corresponding clusters 
are more anisotropic on the square lattice than on the triangular lattice. However, both the 
hulls contain chiral dangling ends and look very similar on small length scales. 



III. SCALING RELATIONS 



Analogous to the cluster related quantities, the hull related quantities also exhibit critical 
properties at the percolation threshold Pc- The corresponding critical exponents associated 
to hull related quantities are expected to satisfy certain scaling relations among themselves 
similar to the scaling relations satisfied by the critical exponents of the cluster related quan- 
tities. The scaling theory of the cluster related quantities for the DSP clusters has already 
been developed in Ref . ^ . In this section, the critical exponents of the analogous hull related 
quantities will be defined and their scaling relations will be developed. The size of a hull H 
is given by the number of occupied sites present on the hull as the size of a cluster S is given 
by the number of occupied sites in a cluster. At Pc, the size of the infinite cluster goes as 

^oo ~ L'^f (1) 



with the system size L and df is the fractal dimension of the infinite c 



uster. The value of 

df has already been estimated on both the square 1] and triangular {2] lattices. For hulls, 
a similar relation is assumed here between the size of the largest hull with the system 
size L as 

~ L''" (2) 
where dn is the hull dimension. The value of dn then can be determined measuring H, 



for different system sizes L. It may be noted here that the so-called Ziff's method 



IJ is 



not suitable to determine the hull dimension of the anisotropic clusters generated here. In 



Ziff's method, the hnear distance is measured for a given number of sites along the hull. 
Since there are two length scales (.^n and involved in the DSP clusters, there will be a 
crossover from C,± to ^|| as one measures the linear distance from smaller number of sites to 
larger number of sites. The average length of a given number of occupied sites will be then 
strongly dependent on the initial point chosen. 

From equations ^ and |2l one can easily verify the relation between the cluster size Sqq 
and the corresponding hull size Hoo for the spanning clusters and it is given by 

~ x = dn/df (3) 

where dn and df are the fractal dimensions of the hull and the corresponding percolation 
cluster, respectively. The value of dn obtained from equation El then could be verified 
measuring x and df independently. 

The hull size distribution is defined as 

Ph{p) = Nn/Ntot (4) 

where Nh is the number of hulls of size H and Ntot is the total number of hulls generated same 
as the number of clusters generated. Analogous to the form of the cluster size distribution 
function, the scaling function form of the hull size distribution is assumed, as 

Ph{p) = H--"+^[H"^{p-p,)] (5) 

where th and (Th are two exponents. The same scaling function form has already been 
verified for the spiral percolation hull[lo|. 

Different moments of the hull size distribution Ph{p), Y^'h H^Ph{p) are expected to be 
singular as p Pc- The primed sum represents the sum of all finite hulls. The first, 
second and third moments xh, x'h ^^"^ Xh '^^ Ph{p) are calculated. The first moment 
Xh = J2'h HPh{p) is the average hull size. The moments xh, x'h Xh diverge with their 
respective critical exponents •jh, and rjH ai p = Pc- The critical exponents are defined as 

Xh ^ \p - Pc\~^" , x'h ^ \P - Pc\~^" ^ Xh ^ \P - Pcl'"^" ■ (6) 

Since the hull related quantities are just different moments of the hull size distribution 
function Pnip), the critical exponents 'Jh, ^h and rjn then should be related to th and 



an, exponents related to Pnip)- It can be shown that the kth moment of the hull size 
distribution become singular as 

The following scaling relations then can easily be obtained putting appropriate values of 
k, order of the moment of Ph{p), in equation [3 and one obtains, 

7h = (3 - th)/<jh, Sh = {4:- th)/<jh t r]H = {5 - th)/<jh- (8) 

Eliminating th and an from equation |H1 a scaling relation can be obtained as 

rjH = 25 H - IH- (9) 

The exponents th and au can also be estimated from the measured exponents 7^, 5h and 
rjH using the following relations 

Th = (35// - A'^h)/{5h - 1h) = {'^Vh - 56H)/{r]H - Sh) = {^Vh - ^1h)/{vh - 1h), ^^^^ 
an = 1/ {Sh - 1h) = 1/ {vh - Sh) = 2/ {t]h - 7h)- 

The values of the hull fractal dimension and other critical exponents will be estimated 
below and the scaling theory will be verified. 



IV. RESULTS AND DISCUSSIONS 



The critical properties of the DSP clusters have already been studied through Monte 
Carlo simulation and finite size scaling on both the square and triangular lattices P,!^]- The 
values of the fractal dimension and other critical exponents of the cluster related quantities 
were estimated numerically and are listed in table HI The values of the critical exponents 
were found different from that of the other percolation models like OP, DP and SP and 
consequently the DSP model belongs to a new universality class[l]. It was observed that 
the spanning clusters are less anisotropic and more compact on the triangular lattice than 
those on the square lattice. This has happened due to the extra flexibility in the spiraling 
constraint on the triangular lattice. The change in the shape and compactness of the clusters 
had significant effect on the values of some of the critical exponents and fractal dimension. 
They were found considerably different on these two lattices for the DSP model and leads 



to a breakdown of universality between the square and triangular lattices in 2D [2]. Below, 
estimates of the hull fractal dimension and other critical exponents of the hull related quan- 
tities will be presented for anisotropic DSP cluster hulls for the first time. A comparison 
will be made between the results obtained on the square and triangular lattices, and also 
with that of the other percolation models. 

To extract hulls, clusters are generated both on the square and triangular lattices at their 
percolation thresholds pc ~ 0.655 and 0.570 respectively. Different lattice sizes starting from 
L = 2^ to 2^^ have been used. Results are obtained through both single lattice (2^^ x 2^^) 
Monte Carlo and finite size scaling methods. A total number Ntot of 5 x 10'^ spanning 
clusters are generated for each lattice size L and every p, site occupation prbability. Since 
each cluster has an associated hull, the number of hulls extracted are also Ntot- The size of 
the hull H is given by the number of occupied sites belonging to the hull. 

The hull fractal dimension of the DSP clusters is determined following the scaling relation. 
Hoc ~ L'^" , equation 121 In Figure IH average hull size of the spanning clusters Hoc is plotted 
against the system size L. The squares represent the square lattice data and the triangles 
represent the triangular lattice data. The values of the hull dimension dn are obtained as 
dn = 1.458 ± 0.008 for the square lattice and dn = 1.463 ± 0.004 for the triangular lattice. 
The hull fractal dimension dn has also been measured by box counting method generating 
spanning clusters on the largest lattice (2^^ x 2^^) considered here. The number of boxes 
occupied with a hull site Nsie) ~ e*^^, where e is the box size. The results are shown in 
the inset of Figure lU The values of dn obtained in the box counting method are given as 
1.44 ± 0.01 for the square and 1.45 ± 0.01 for the triangular lattice. The errors quoted for 
both the methods are the least square fit errors taking into account the statistical error of 
each data point. The results obtained through box counting method and finite size scaling 
are within error bars. There are few things to notice. First, the hull fractal dimensions 
dn measured on the square and triangular lattices in 2D are found the same within error 
bars, approximately 1.46, whereas the cluster fractal dimension df were found different on 
the same lattices|2]. It could also be seen that, the hulls shown in Figure |21 and El of the 
clusters generated on the square and triangular lattices respectively are very identical on 
small length scales. However, DSP clusters are more compact and less anisotropic on the 
triangular lattice than on the square lattice. As a consequence, the cluster fractal dimension 
df and other critical exponents of the DSP clusters were found different on these two lattices 



(see Table H}. Thus, the extra flexibihty given on the rotational constraint on the triangular 
lattice was only able to modify the critical behaviour of the whole cluster but unable to 
modify the critical properties of the external perimeter, the hull. 

Second, the hull fractal dimension dn obtained here is smaller than that of OP cluster 
hulls (7/4). The OP cluster hulls generally contains long fjords or bays which are absent 
here. For OP cluster hulls, the value of dn = 7/4 was conjectured by Sapoval a/ js] 
through a relation dn = 1 + l/iy connecting the hull dimension dn and the connectivity 
length exponent u studying the diffusion fronts. (In the case of OP, the value of u is 
4/3). This prediction was supported through large scale simulation by Ziff and also 
proved analytically by Saleur and Duplantier jl^. It has already been observed that the 
conjecture does not hold true in the case of SP cluster hulls. The values of the connectivity 
exponents and hull fractal dimensions obtained for the SP clusters are: v{SP) ~ 1.116 
and dui^SP) ~ 1.476 {1 + \/u ~ 1.896) on the square lattice and i^{SP) ~ 1.136 and 
dniSP) ^ 1.466 (1 + 1/u ^ 1.880) on the triangular lattice Moreover, the DSP 

clusters are anisotropic and there are two connectivity length exponents z/y and z/_l. The hull 
dimension dn here then may not be related to the connectivity length exponents z/y and z/_|_ 
of the anisotropic clusters in a simple manner as it was predicted by Sapoval et al for the 
case of isotropic clusters generated in OP. Interestingly, it is found that the hull dimension 
dn calculated form the following relation 

dH = l + (11) 



is very close to the numerical values obtained here. Using equation [TTJ the hull dimensions 
dn are obtained as dn ~ 1-46 on the square lattice and du ~ 1.47 on the triangular lattice. It 
seems that for the anisotropic clusters generated in DSP, the hull dimension dn is connected 
to the connectivity length exponents v\\ and by the proposed relation in equation ITT] and 
the relation proposed by Sapoval et al is valid only for the isotropic clusters generated in 
OP. However, an argument similar to the one given by Sapoval et a/[8] would be possible 
only from the study of "anisotropic diffusion fronts" . 

Third, the values of the hull fractal dimensions obtained for both the DSP and SP clusters 
are close to the fractal dimension = 4/3 of the externally accessible perimeter defined 
by Grossman and AharonylS]. The external perimeter defined by Grossman and Aharony 
includes the sites available to a finite-size particle, coming from the outside, that is touching 



the occupied sites on the cluster. This external perimeter excludes deep fjords or bays and 
consequently the fractal dimension found 4/3. However, in the case of SP and DSP clusters, 
the hull is compact (or smooth) in comparison to the OP cluster hull due to the presence 
of rotational constraint in these models. The rotational constraint produces compact chiral 
dangling ends on the external boundary and makes the hull fractal dimension dn closer 
to that of the Grossman-Aharony external perimeter. It should be mentioned here that 
similar value of hull fractal dimension close to 4/3 was also obtained numerically by Meakin 
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percolation fronts by Birovljev et a/[l^ and in the corrosion of thin Aluminum film by 
Balazsjiol. Deviation from the fractal dimension of 7/4 was also found in the study of self- 
stabihzed etching of random systems by Sapoval et a/|2^. Recently, Aizenman et a/ proved 
exactly the fractal dimension of Grossman-Aharony external perimeter is 4/3 21 1. 

The relation between the hull dimension du and the fractal dimension dj oi the per- 
colation clusters is given by equation 01 Hoo ~ S^,x = dn/dj. This is verified here and 
the exponent x is calculated. In Figure the average hull size Hoo is plotted against the 
corresponding average size of the spanning clusters 6*00 for both the square and triangular 
lattices. The slope x = du/dj is found as x = 0.839 ± 0.006 on the square lattice and 
X = 0.820 ± 0.006 on the triangular lattice. The errors are least square fit error taking 
into account the statistical error of each data point. The values of x are found close but 
only slightly different on the two lattices. This small difference in the value of x is consis- 
tent with the small difference in the cluster fractal dimensions df on the two lattices. The 
hull dimension dn could be estimated from the relation dH_= x x df and compared with 
the measured values. Taking the values of df from Ref. O], given in Table HJ it is found 
that dn = 1-45 ± 0.01 on the square lattice and dn = 1-46 ± 0.01 on the triangular lattice 
which are within error bars of the measured value ~ 1.46. The errors quoted here are the 
propagation errors. 

Now, the critical properties of different moments of the hull size distribution function are 
studied. The critical exponents related to the different moments of the hull size distribution 
Ph{p) are already defined in section III and the scaling relations are also described there. In 
this section, the values of the critical exponents are estimated and the scaling relations are 
verified. The first three moments xh^ x'h ^"^^ x"h plotted against \p — pc\ for the square 
lattice in Figure IHl^a) and for the triangular lattice in Figure Effe). The system size is taken 



as L = 2048. The circles represent the average hull size xh, the squares represent x^, the 
second moment and the triangles represent x'hi the third moment in both the plots. The 
values of the exponents obtained are 7^ = 1.68±0.02, 5h = 3.66±0.03 and rjH = 5.70±0.05 
for the square lattice. For the triangular lattice, the values of the exponents obtained are 

= 1.71 ± 0.02, 5h = 3.69 ± 0.03 and rju = 5.73 ± 0.05. The errors quoted here are the 
standard least square fit error taking into account the statistical error of each single data 
point. The system size L dependence of the values of the critical exponents of different 
moments of the hull size distribution on the square and triangular lattices is also studied 
measuring the exponents' values at different L. The values of •jh, and rjn are plotted 
against the inverse system size 1/L in Figure [7|^a) , (6) and (c) respectively. The squares 
represent the square lattice data and the triangles represent the triangular lattice data in 
each plot. The values of the exponents obtained on these two lattices are converging to the 
values of the exponents obtained for L = 2048 through MC simulation. The values of the 
hull moment exponents are also within the error bars on the square and triangular lattices. 
The hull fractal dimension dn is already found the same on these two lattices. It should 
be mentioned here that the values of the critical exponents of the analogous cluster related 
quantities were found different on the square and triangular lattices Q]. Thus, breakdown 
of universality occurs in the cluster properties whereas universality holds for the associated 
hull properties in the DSP model in 2D between the square and triangular lattices. This is a 
new observation and appears for the first time in a percolation model. The scaling relations 
could be verified now. The value of 26h — Ih = 5.64 ± 0.08 is very close to the value of the 
exponent rjn ~ 5.70 for the square lattice. For the triangular lattice, 26h — 1h = 5.67 ±0.08 
and rjH ~ 5.73. The scaling relation rjH = 26h — 1h (equation IH)) is then valid within the 
error bars on both the lattices. 

In table m the values of the hull fractal dimension and other critical exponents are 
summarized and a comparison of the results obtained on the square and triangular lattices 
as well as on the SP model has been made. It can be seen that the universality holds for the 
hull results on the square and triangular lattice, as it is already mentioned, though there is 
a breakdown of universality for the cluster properties on the same lattices. It is interesting 
to note that the values of the exponents obtained here, especially the hull fractal dimension 
dn, are very close to that of the SP (percolation under rotational constraint only) clusters 
hulljl(l| on both the lattices. However, the SP clusters were found much more compact than 



the DSP clusters and are isotropic. The fractal dimension ci/ of SP clusters on the square and 
triangular lattices are ~ 1.957 and 1.969 respectively[lQ:] whereas that of the DSP clusters 
are ~ 1.733 [Ij and 1.775 [3] respectively. It can be seen from the hull structures, given in 
Figure El and El that there exist chiral (clockwisely rotated) dangling ends on the perimeter 
as in the case of SP cluster hulls. The hull properties of both the DSP and SP clusters are 
then determined by the existence of chiral dangling ends on the external perimeter. These 
chiral dangling ends are generated due to the presence of the rotational constraint in both 
the models. The directional constraint then has a little effect on the hull critical properties 
and it is mostly determined by the rotational constraint present in the model. In a physical 
situation like transport of classical charged particles in disordered systems in the presence of 
crossed electric and magnetic fields, the magnetic properties then could be extracted from 
the external perimeter only whereas the electrical properties could be obtained from the 
bulk of the cluster. 

Finally, the form of the hull size scaling function Ph{p) = H^'^'^^^f[H'^^ {p—pc)] is verified. 
The values of the exponents th and cxh have been estimated using the scaling relation 
given in equation ^1 For the square lattice, the estimates of th and an are obtained 
as Th = 2.17 ± 0.02 and an = 0.498 ± 0.003. These values are within the error bars 
of the values of th and an obtained on the triangular lattice as th = 2.16 ± 0.02 and 
an = 0.497 ± 0.003. The errors mentioned correspond to propagation error only. The 
scaling function form is verified through data collapse by plotting Ph{p)/Ph{Pc) against the 
scaled variable H'^" {p — pc) for both the square and triangular lattices in Figure IHfa) and 
(6), respectively. The hull size varies from 64 to 16384 and {p — Pc) is in the range 0.007 
to —0.06. Not only a reasonable data collapse is observed for both the lattices but also the 
scaling function forms are found identical on the two lattices. 

V. CONCLUSION 

Hulls of the DSP clusters are extracted generating large clusters at the percolation thresh- 
old and their properties are analyzed on both the square and triangular lattices in 2D. Re- 
sults are compared with that of other percolation models. The values of the hull fractal 
dimension dn and critical exponents •jh, and rjn related to the hull size distribution 
function Ph{p) are determined. It is found that the values obtained for the hull fractal di- 



mension and the critical exponents related to the hull size distribution are the same within 
error bars on the square and triangular lattices unlike the cluster fractal dimension df and 
cluster related critical exponents of the DSP model. Thus, the universality of the hull critical 
exponents holds for the square and triangular lattices in 2D whereas there is a breakdown 
of universality in the corresponding cluster properties. The value of dn ~ 1-46 obtained 
here is close to 4/3, the fractal dimension of the Grossman- Aharony external perimeter and 
away from 7/4 proposed by Sapoval et a/ for the ordinary percolation hull. A new conjecture 
is made relating dn and two connectivity length exponents z/y and u± for the anisotropic 
clusters. It is found that the values of dn and other hull critical exponents are close to that 
of the SP cluster hulls. Hull critical properties are then mostly determined by the rotational 
constraint than the directional constraint present in the model. The critical exponents esti- 
mated here satisfy the scaling relations among themselves within error bars. The assumed 
scaling function form has been verified and found similar on both the square and triangular 
lattices. 
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TABLE L Numerical estimates of the critical exponents and fractal dimension measured for the 
DSP clusters on the square and triangular lattices. Some of the critical exponents and the fractal 
dimension are significantly different on the two lattices. 



Lattice Type 


Model 


du 


d% 


X 


dn = xdf 


IH 


5h 


r]H 


Square 


DSP 


1.458 


1.46 


0.839 


1.45 


1.68 


3.66 


5.70 






±0.008 


±0.02 


±0.007 


±0.01 


±0.02 


±0.03 


±0.05 






dH = 1 


44 ±0.01 


(Box Counting) 










SP 


1.476 




0.74 




1.82 


3.75 








±0.005 




±0.02 




±0.01 


±0.02 




Triangular 


DSP 


1.463 


1.47 


0.820 


1.46 


1.71 


3.69 


5.73 






±0.004 


±0.02 


±0.006 


±0.01 


±0.02 


±0.03 


±0.05 






du = 1 


45 ±0.01 


(Box Counting) 










SP 


1.466 




0.76 




1.91 


3.87 








±0.016 




±0.02 




±0.01 


±0.03 





TABLE II: Comparison of the hull exponents and hull fractal dimension measured for the DSP 
model on the square and triangular lattices as well as with that of the SP model, is the hull 
fractal dimension obtained through the proposed conjecture: dn = 1 ± i^±/{i^\\ ± i^±)- dn is also 
measured from the relation du = xdf considering the values oi df as 1.733 and 1.775 (table I) 
for the square and triangular lattices respectively. Measurements are found consistent in different 
methods adopted. The exponents are almost same within the error bars on the two lattices. The 
hull fractal dimension in the DSP model is found within error bar of that of the SP model. 
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FIG. 1: Selection of empty nearest neighbours for occupation in a MC time step is demonstrated 
here on the square (a) and triangular (b) lattices. Black circles are the occupied sites and the 
open circles are the empty sites at a given time. Thick long arrows from left to right represent the 
directional constraint. The presence of clockwise rotational constraint is shown by encircled dots. 
The central site is occupied from site 2 and shown by short thick arrows. On the square lattice, 
site 3 on the right is always eligible for occupation due to the directional constraint, marked by 
dotted arrow and sites 4 and 1 are eligible for occupation due to the rotational constraint, marked 
by thin arrows. Similarly on the triangular lattice, site 4 due to the directional constraint and sites 
5, 6 and 1 due to the rotational constraint are eligible for occupation. 



FIG. 2: Hull of a spanning cluster of size S = 23964 on a 256 x 256 square lattice at the percolation 
threshold Pc = 0.655. The size of the hull is H = 7938. It could be seen that the dangling ends are 
clockwisely rotated and the cluster is anisotropic. 



FIG. 3: Hull of a spanning cluster of size S = 15648 on a 256 x 256 triangular lattice at the 
percolation threshold Pc = 0.570. The size of the hull is H = 4370. There are chiral dangling ends 
but the cluster is less anisotropic than the cluster on the square lattice. 




FIG. 4: Average hull size H^o of the spanning clusters versus the system size L. The squares 
represent the square lattice data and the triangles represent the triangular lattice data. The hull 
fractal dimensions are obtained as dn = 1.458 ± 0.008 for the square lattice and 1.463 ± 0.004 for 
the triangular lattice. In the inset, the number of boxes NB{e) is plotted against the box size e. 
The hull fractal dimensions obtained by box counting method are dn = 1.44 it 0.01 for the square 
lattice and 1.45 it 0.01 for the triangular lattice. The values obtained in two different methods are 
within error bars. The value of dn is also within error bars on the two lattices. 
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FIG. 5: Plot of average hull size versus corresponding average cluster size Soo of the spanning 
clusters. The squares represent the square lattice data and the triangles represent the triangular 
lattice data. The values of the exponent x are obtained as x = 0.839 ± 0.007 for the square lattice 
and X = 0.820 it 0.006 for the triangular lattice. 
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FIG. 6: Plot of the first, second and third moments xh, x'h Xh ^^^^ ^i^^ distribution versus 
|p — Pel for the lattice size L = 2048 on the square (a) and triangular (6) lattices. Different symbols 
are circles for xh, squares for x'h triangles for x'h each lattice. The solid lines represent 
the best fitted straight lines through the data points. The corresponding critical exponents are 
found as jh = 1-68 ± 0.02, 5h = 3.66 it 0.03 and rjH = 5.70 =b 0.05 for the square lattice and 
7if = 1.71 ± 0.02, Sh = 3.69 ± 0.03 and rjH = 5.73 ± 0.05 for the triangular lattice. 
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FIG. 7: Plot of the hull moments exponents (o), 5h (6)and ■qH (c) against the inverse system 
size 1/L. The system sizes considered are: L = 128, 256, 512, 1024 and 2048. The squares represent 
square lattice data and the triangles represent triangular lattice data. As L — > oo, the exponents 
are converging to the values corresponding to L = 2048. The values of the exponents are within 
the error bar on the two lattices. 




h^^Cp-p,) 

FIG. 8: Plot of scaled hull size distribution Ph{p)/ Ph{Pc) versus scaled variable H'^"{p — Pc) 
for different values of p on the square (a) and triangular (b) lattices. The value of ajj is taken 
as an = 0.498 (square) and an = 0.497 (triangular). The hull size H changes from 64 to 16384. 
The data plotted correspond to p - Pc = 0.007(x), 0.005(0), -0.035(V), -0.04(0), -0.045(<), 
— 0.05(*), —0.055(A), —0.06(0) for both the plots. Reasonable data collapse are observed on both 
the lattices. Since the value an is almost the same, the scaling function form looks almost identical 



